/* Dorussen and Ward Data (+MTOPS) */
/* This file replicates Figures 1,3,5&6 in the paper */
/*You will need to change the directories below to account for your directory structure*/


/* model #1 from final paper (pg 45)*/
%INCLUDE "H:\My Documents\My SAS Files\9\Macro_TVEM_v210\Macro_TVEM_v210_logistic.sas";
libname work "H:\My Documents\My SAS Files\9\"; 

%TVEM_logistic (
  mydata = work.fulldata_fatalmid_small, 
  id = dyadid, 
  time = year, 
   dep = force, 
   tcov = intercept i2_numigo, 
   cov_knots = 6 6,  
   cov = low_polity2 noncontiguity low_depend_ksg2,
   evenly = 0,  
   scale = 100,  
  oddsratio = 0, 
  outfilename = model_1_finaldraft_p45_fatalmid.csv);      
 
 /***************************************/

/*model #3 from the final paper */
%INCLUDE "H:\My Documents\My SAS Files\9\Macro_TVEM_v210\Macro_TVEM_v210_logistic.sas";
libname work "H:\My Documents\My SAS Files\9\"; 

%TVEM_logistic (
  mydata = work.fulldata_fatalmid_small, 
  id = dyadid, 
  time = year, 
   dep = fatalmid, 
   tcov = intercept i2_numigo, 
   cov_knots = 6 6,  
   cov = low_polity2 noncontiguity low_depend_ksg2,
   evenly = 0,  
   scale = 100,  
  oddsratio = 0, 
  outfilename = model_3_finaldraft_p47_fatalmid.csv);    
  
/***************************************/

/* model #5 from final paper */
%INCLUDE "H:\My Documents\My SAS Files\9\Macro_TVEM_v210\Macro_TVEM_v210_logistic.sas";
libname work "H:\My Documents\My SAS Files\9\"; 

%TVEM_logistic (
  mydata = work.fulldata_fatalmid_small, 
  id = dyadid, 
  time = year, 
   dep = force, 
   tcov = intercept i2_num_mtops, 
   cov_knots = 6 6,  
   cov = low_polity2 noncontiguity low_depend_ksg2,
   evenly = 0,  
   scale = 100,  
  oddsratio = 0, 
  outfilename = model_5_finaldraft_p49_fatalmid.csv);    
    
/***************************************/

/* model #6 from final paper */
%INCLUDE "H:\My Documents\My SAS Files\9\Macro_TVEM_v210\Macro_TVEM_v210_logistic.sas";
libname work "H:\My Documents\My SAS Files\9\"; 

%TVEM_logistic (
  mydata = work.fulldata_fatalmid_small, 
  id = dyadid, 
  time = year, 
   dep = fatalmid, 
   tcov = intercept i2_num_mtops, 
   cov_knots = 6 6,  
   cov = low_polity2 noncontiguity low_depend_ksg2,
   evenly = 0,  
   scale = 100,  
  oddsratio = 0, 
  outfilename = model_6_finaldraft_p50_fatalmid.csv);    
    
/***************************************/
